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\ The quantum-mechanical framework in which observables are associated 

with Hermitian operators is too narrow to discuss measurements of such im- 
portant physical quantities as elapsed time or harmonic-oscillator phase. We in- 
i-£h I troduce a broader framework that allows us to derive quantum-mechanical limits 

on the precision to which a parameter — e.g., elapsed time — may be determined 
via arbitrary data analysis of arbitrary measurements on N identically prepared 
quantum systems. The limits are expressed as generalized Mandelstam-Tamm 
uncertainty relations, which involve the operator that gener ates displacements of 
the parameter — e.g., the Hamiltonian operator in the case of elapsed time. This 
approach avoids entirely the problem of associating a Hermitian operator with 
the parameter. We illustrate the general formalism, first, with nonrelativistic un- 
■ certainty relations for spatial displacement and momentum, harmonic-oscillator 

phase and number of quanta, and time and energy and, second, with Lorentz- 
invariant uncertainty relations involving the displacement and Lorentz-rotation 
parameters of the Poincaré group. 
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1 Introduction 



The goal of quantitative experiments in physics is to determine a set of parameters 
to some level of confidence. In general this determination entails complex methods 
of data analysis applied to observed data. From this point of view, the conventional 
description of measurements in quantum theory, tied to the use of Hermitian opera- 
tors to represent observable quantities, provides too narrow a framework, because for 
many experimental parameters, time being an example, there is no suitable Hermitian 
operator. 

In this paper we employ a broader framework for describing the quantum-mechan- 
ical determination of parameters such as time M. In this framework measurements 
are described in the most general way permitted by quantum mechanics — in terms 
of so-called "positive-operator-valued measures" (POVMs). The role of a quantum 
measurement is to provide data from which one infers the parameter of interest by 
classical methods of parameter estimation. There is no need to associate a Hermitian 
operator with the parameter, and generally there is no such Hermitian operator. We 
derive quantum restrictions on determining a parameter by considering optimal mea- 
surements and optimal methods of parameter estimation. The quantum restrictions 
are stated as uncertainty relations that involve the parameter and the operator that 
generates displacements of the parameter, time and the Hamiltonian operator being 
an example. 

Uncertainty relations are central to the interpretation of quantum theory, yet in 
many cases of interest it is impossible to formulate an uncertainty relation if one insists 
that both quantities have an associated Hermitian operator. Hilgevoord and Uífink 
give an excellent summary of the defects of Standard uncertainty relations and of 
the motivation for parameter-based uncertainty relations. Mandelstam and Tamm || 
derived the first parameter-based uncertainty relation, for time and energy, by treating 
elapsed time as a parameter to be determined by measurement of a conventional 
observable that varies with time. Helstrom [H] and Holevo [H] pioneered the modern 
study of parameter-based uncertainty relations, by considering quantum restrictions 
on how well one can determine a parameter from the results of general quantum 
measurements described by POVMs. Other authors [p], 0, || |7|, [8j have formulated 
parameter-based uncertainty relations in various contexts. 

Here we present a general theory of parameter-based uncertainty relations and 
explore in some detail the question of finding optimal quantum measurements that 
achieve the lower bound set by the uncertainty relation. We devote Section ^| to 
summarizing the framework for quantum parameter estimation and the corresponding 
generalized parameter-based uncertainty relations. Section [O] develops the general 
theory for mixed quantum states (density operators). Section |2.2| specializes the gen- 
eral theory to pure states that are generated by a single-parameter unitary operator, 
a case that occupies the remainder of the paper. In Section |3| we develop a general 
description of global optimal measurements that saturate the lower bound in the gen- 
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eralized uncertainty relation. Section |] illustrates the parameter-based uncertainty 
relations with various examples of nonrelativistic uncertainty relations: spatial dis- 
placement and momentum in Section |4.1| , harmonic-oscillator phase and number of 
quanta in Section p~2| , and time and energy in Section p~3j Section || applies the 
parameter-based uncertainty relations to the displacement and Lorentz-rotation pa- 
rameters of the Poincaré group, leading ultimately to relativistically invariant uncer- 
tainty relations for the invariant space-time interval of special relativity and the boost 
and spatial-rotation parameters of Lorentz transformations. Section ^| concludes with 
a brief discussion. 



2 Generalized Uncertainty Relations 
2.1 Uncertainty relations for mixed states 

Consider N replicas of a quantum system. Each replica is prepared in the same 
quantum state (density operator) p(X), which is parametrized by the single parameter 
X. In the following a subscript X on an expectation value denotes an expectation 
value with respect to p{X). Braunstein and Caves ÜJ consider a general smooth path 
on the space of density operators, 

p(x) = J2 Pj \j)(j\, (i) 

j 

where both the eigenvalues pj and the eigenvectors \j) can change along the path. A 
path is specified by giving the tangent vector 

The Hermitian operator h, which can depend on X, generates the infinitesimal changes 
in the eigenvectors of p(X): 

p(X + dX) = P (X) + p\X)dX = J2(P 3 + d Pj )e- idxh \j)(j\e idxh . (3) 

i 

Notice that h can be replaced by 

Ah = h- (h) x (4) 

in Eqs. (0) and (§) without changing the path. 

The most general measurement permitted by quantum mechanics [KJ can be 
described by a set of bounded, non-negative, Hermitian operators E(Ç)d£ (generaliza- 
tions of projection operators), which are complete in the sense that 

J d^ È(Ç) = í = (unit operator) . (5) 
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The quantity £ labels the "results" of the measurement; written here as a single 
continuous real variable, it can be discrete or multivariate. The operators E(^)d^ 
make up what is called a "positive-operator-valued measure" (POVM). The probability 
distribution for result £, given the parameter X, is 

p(t\X) = tr(È(Op(X)) . (6) 

The properties of the POVM are just those needed to make p(Ç\X) a normalized 
probability distribution. 

Let £1, . . • , £jv denote the results of measurements on the N replicas of our quantum 
system. A general form of data analysis uses a function 

X cst = X cst (£i, . . . , £at) (7) 

to generate an estimate X est for the parameter X, based on the data £1, . . . ,£jv ob- 
tained from the measurements and nothing else. 

To characterize how precisely the N measurements are able to determine the pa- 
rameter X, we need something a bit more complicated than the obvious choice, the 
variance of the estimator, ((AX est ) 2 )x = ((X est — (X est )x) 2 )x- The reason is that 
the variance does not take into account two important possibilities. First, even if the 
estimator has a small variance, it might be systematically biased away from the true 
parameter value — i.e., (X est ) x might not equal X — and thus give a poor estimate. 
Second, the estimator might have different "units" from the parameter, thus making 
it difficult to intèrpret the variance of the estimator as a measure of precision in de- 
termining X. Both the amount of bias and the difference in units can depend on the 
parameter, i.e., on location along the path. To remedy these difficulties, we quantify 
the estimate's deviation from the parameter by |J 

ÓX = ■ - Xcst ; - X . (8) 

\d(X cst ) x /dX\ y ' 

The derivative d(X est )x/dX removes the local difference in the "units" of the esti- 
mator and the parameter, and then the units-corrected estimator is compared to the 
parameter X, not to the mean value of the estimator. As a statistical measure of the 
precision of the estimation, we use the second moment of SX. 
There is a lower bound on the second moment of SX: 

Braunstein and Caves [0] derive the ultimate lower bound in two steps, in contrast 
to derivations by Helstrom |j| and Holevo [5(Chap. VI. 2)], both of whom proceed to 
the ultimate lower bound in a single step that obscures the conditions for achieving 
the ultimate lower bound. The two steps in the Braunstein-Caves derivation are 
displayed as the two inequalities in Eq. (P). The first inequality is a bound that 
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applies to all estimators X est for a fixed probability distribution p(Ç\X), i.e., for a 
fixed quantum measurement. The second inequality is a bound that applies to all 
quantum measurements. 

In the first inequality in Eq. @, 



is the Fisher information |TTJ associated with the probability distribution p(Ç\X). 



The first inequality is an expression of the Cramér-Rao bound of classical estimation 
theory ||11|| , which places a lower bound on the variance of any estimator X est that 
is applied to data drawn from the distribution p(Ç\X). An estimator that saturates 
the first inequality in Eq. (||) — and, hence, attains the Cramér-Rao bound — is called 
an efficient estimator. The lower bound in the first inequality can always be achieved 
asymptotically for large N by using maximum-likelihood estimation |T3], but except 
for special distributions, there is no efficient estimator for finite vàlues of N. 

The second inequality in Eq. (|9|) holds for any POVM E(Ç)dÇ. The second in- 
equality is written in terms of a line element ds 2 , which defines a "statistical distance" 
|T3| that measures the distinguishability of neighboring quantum states and provides 



a natural Riemannian geometry on the space of density operators. The explicit form 
that Braunstein and Caves Jï]] (see also [4,5(Chap. VI. 2)]) find for the line element is 



ds 2 /dX 2 = ((Cé(p')) ) x = tr(> X p (p')) , (11) 
where C p is a super-operator that, in the basis that diagonalizes p, takes the form 

£fi(Ó)= £ —^—0 3k \j)(k\. (12) 

If p has no zero eigenvalues, C p is the inverse of the super-operator defined by 1Z P (Ò) = 
\{pÒ + Óp); hence, C p was denoted 1ZJ 1 by Braunstein and Caves [0]. Helstrom ||] 
and Holevo [5(Chap. VI. 2)] call Cp(p') the "symmetric logarithmic derivative" of p. 
The lower bound in the second inequality can be achieved by using a POVM such 
that the operators E(Ç) are one-dimensional projection operators onto orthonormal 
eigenstates of the Hermitian operator Cp(p') 0. The conditions given here and above 
for achieving the two lower bounds in Eq. (|9|) are sufficient, though they are not always 
necessary. 

The line element (|ïï|) on the space of density operators arises here from quantifying 
precisely the quantum restrictions on determining a parameter — and, hence, the quan- 
tum restrictions on the statistical distinguishability of neighboring density operators 
p{X) and p(X + dX). The same line element can also be gotten by defining a natural 
mètric on density operators in terms of the correlation between pairs of conventional 
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observables. The reader interested in this alternative route to the mètric can find it 
spelled out in @, Q, together with references to related work. 

Although not pointed out by Braunstein and Caves, the lower bound (|9|) does not 
improve if one allows measurements that do not factor into separate measurements 
on each of the N replicas. One can see this by treating the N replicas as a single 
composite system with density operator 

pW(x) = p(X)®···®p(X) . (13) 

Applied to this composite system, the bound (§) takes a form {(5X) 2 )x > (dX/ds^) 2 
that holds for all quantum measurements on the product space of the composite sys- 
tem. It is not difficult to show, however, that for iV-replica product states, the line 
element on the product space reduces to N times the single-replica line element, i.e., 

~ix^ = 'iir^iir) ) = Nt <^) • < 14 > 

thus giving the same lower bound as in Eq. @ . This result provides a limited answer 
to a question raised by Peres and Wootters ||15|| : when a composite system is made 
up of replicas all prepared in the same quantum state, can measurements on the 
composite system better distinguish states than can separate measurements on each 
of the replicas? For the very special case of two neighboring states, the answer is no. 

We pause at this point to take stock of what has already been presented. The 
bound @, together with Eq. (|ïï|), is a general species of uncertainty relation, which 
restricts one's ability to determine a parameter from the results of quantum mea- 
surements. This uncertainty relation applies to mixed quantum states, allows for 
measurements that are not described by projection operators, and includes the possi- 
bility of múltiple measurements. On the other hand, precisely because this uncertainty 
relation is so general, we find it instructive in what follows to specialize in three ways, 
thus permitting us to make closer contact with Standard uncertainty relations. 

For the first specialization we assume that the eigenvalues of the density operator 
do not change along the path — i.e., dpj = in Eq. @ — which simplifies Eq. (§) to 

p' = -i[h,p] = -i[Ah,p]. (15) 

This first specialization means that the path is generated by a unitary transformation; 
keep in mind, however, that we still allow the local generator of the transformation, 
h, to depend on X. As a consequence of this first specialization, we can write 

C0) = 2i £ PL·L^Ah jk \j)(k\ = 2rh, (16) 
{JMPj+Pk^o} Pj + P k 

where we introduce 5h as a shorthand for ^C^(p'), and the line element (|TTJ) becomes 

^- = mhf) x = 2£fe+p fe )í^f YWi 2 < M(AhY) x . (i?) 

auX j,k \Pj ~r Pk/ 
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Notice that in this line element we can drop the restriction on the sum, since under 
any procedure for approaching the boundary on which one or more eigenvalues of p 
vanishes, the terms for which pj + = do not contribute. 

A consequence of the last inequality in Eq. (|Ï7| ) is a parameter-based uncertainty 
relation [l,4,5(Chaps. III.2, IV. 7, and VI.3)], 

((5Xf) x ((Ah) 2 ) x > -L , (18) 

which, since it involves the variance of h, resembles Standard uncertainty relations, 
except that X is a parameter and the relation holds for múltiple measurements. The 
corresponding uncertainty relation involving 5h, 

((ÓX) 2 ) x ((fh) 2 ) x >^, (19) 

is stricter [l,5(Chap. VI. 3)], unless equality holds in Eq. ([H]). Equality is equivalent 
to the condition that pjPk\Ahjk\ 2 = for all j and k. In particular, equality holds if p 
is a pure state, but never holds if p has no zero eigenvalues (except in the trivial case 
Ah = 0). 



2.2 Uncertainty relations for pure states 

The second specialization is to assume that p(X) = \4>x){ipx\ is a pure state. This 
assumption implies the first one, which is incorporated in Eq. (|Ï5f) , since a path on 
the pure states must be generated by a unitary transformation. Normalization implies 
that 

. . ( , A -| (ÍM) + (^) fe, - 2 Re(( fe , (%>)) , (20, 

but the freedom to multiply \ipx) by a phase factor means that lm[(ip x \(d\ipx) / dX)} 
can be chosen arbitrarily. Using Eqs. (|Ï5|) and ( p0|) , one can show that 

(^) ± -^-w.xw(^)-^>. (-) 

where (d\ipx) /dX)± is the projection of d\ipx)/dX orthogonal to \ipx)- Equation (|Ï5|) 
can now be written as 

A convenient phase choice, 

Ux\$f)=-i<K)z, (23) 
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leads to a Schròdinger-like equation for \ipx)' 



dtyx) 
dX 



-ih\ip x ) ■ 



(24) 



Notice that the phase freedom in \ipx) is equivalent to the freedom to add a múltiple 
of the unit operator to h. 

Applying our second assumption to Eq. (|Ï6l) , one finds that 



5h 



-i[Ah,p) 



' d\il>xY 
dX 



0>x\ + \ipx) 



x 



dX 



P 



(25) 



Thus, for pure states, the line element (|Ï7|) for statistical distance reduces to 



ds 2 
dX 2 



'd{i>. 



x 



dX 



' d\ipxY 
dX 



A((Ah) 2 ) 



x 



(26) 



which implies, as indicated above, that we can restrict our attention to the uncertainty 
relation (|I|). 

One expects statistical distance, which measures the distinguishability of states, 
to be related to the inner product and thus to the Hilbert-space angle between pure 
states. The square of the infinitesimal Hilbert-space angle d9 between neighboring 

states \if) X ) and \ïp x +dx) is 

dO 2 = [cos-\\^ x \^x +d x)\)} 2 = 1 - \(^x\^x +d x)\ 2 ■ (27) 
The line element d9 2 defines a natural Riemannian mètric, called the Fubini-Study 



mètric [IÜ , |TT| , I5| , on the manifold of Hilbert-space rays. Using Eq. fl2"0|) and the 
further consequence of normalization, 







1 d 2 



2dX- 



;{ipx\^x) 



'd{ip x \ \ ( d\ipxY 



dX 



dX 



+ Re (ipx 



d 2 \^ x ) 
dX 2 



(2í 



one finds that 



de 2 



dX 2 



dX 2 



'd(^ x \ \ ( d\ip x )' 



\ dX 

djipx]* 
dX 



{ dX 

' d\i> x y 

dX 



-ds 2 



' d^xY 
dX 



(29) 



which means that Hilbert-space angle is half the statistical distance defined here. 

The third and final specialization is to assume that the generator h is independent 
of X. This assumption allows us to integrate Eq. ( pip and to write the path on Hilbert 
space as being generated by a single-parameter unitary operator, 



\1>x) 



-iXh 



|-0o> 



(30) 
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where |-0o) is a fiducial state at X = 0. Moreover, this assumption guarantees that 
the expectation value of any function of h is independent of X; thus we can drop 
the subscript X from the mean and variance of h. In particular, we can write the 
uncertainty relat ion (0) as 

((5Xf) x ((Ahf) > -L . (31) 

It is this parameter-based uncertainty relation that occupies us for the remainder of 
this paper. As noted above, this uncertainty relation resembles the Standard uncer- 
tainty relation, except that the relation holds for múltiple measurements and X is a 
parameter, not necessarily corresponding to any Hermitian operator. 



3 Global Optimal Measurements 
3.1 General considerat ions 

The chain of inequalities leading to the uncertainty relation ( |3~ID consists of the two in- 
equalities leading to the statistical distance in Eq. (§) and the inequality involving the 
generator h in Eq. ([T7|). The first inequality in Eq. ([J) can be saturated asymptotically 
for large N by use of maximum-likelihood estimation, and the inequality in Eq. (fl7[) is 
saturated for pure states. Thus the question of achieving equality in the uncertainty 
relation (|3Ï|), provided one allows for many measurements N, reduces to finding an 
optimal measurement, i.e., one that saturates the second inequality in Eq. (H). Notice 
that since the variance of h is constant as a consequence of our third assumption, 
optimal measurements lead to a maximum Fisher information that is constant along 
the path. 

As indicated above, one such optimal measurement uses a POVM such that the 
operators E(£) are one-dimensional projection operators onto orthonormal eigenstates 



of the Hermitian operator Cp(p') = 25h [cf. Eq. (25)]. This measurement has the 
defect, however, that it generally depends on X, thus requiring one to know the value 
of the parameter one is trying to estimate before choosing the optimal measurement. 
Our goal here is to find a global measurement, independent of X, that is optimal all 
along the path. We seek such a global optimal measurement in terms of a POVM 
E(x) dx, where the measurement results are labeled by a single real number x that 



has the same range of vàlues as X. As we discuss further in Section |3.2| , we can hope to 
find an optimal measurement of this form only if the generator h is non-degenerate; if 
the spectrum of h has degenerades, an optimal measurement must acquire information 
beyond that which can be described by a single real number. 

The POVM É(x) dx must, of course, be complete, which means that 

i = / dxÈ{x) . (32) 
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The probability density for result x, given the parameter X, is 



p(x\X) = (ip x \È(x)\ip x ) = (Me iXh È(x)e- iXh \ip ) . (33) 

As noted above, global optimal measurements lead to a Fisher information ( |IÜD that 
is independent of X. This suggests that we require that p(x\X) be a function only of 
x — X, which means the POVM must satisfy a "displacement" property 

é xh È(x)e- ixk = È{x - X) . (34) 

Measurements that satisfy properties ( |3"2| ) and are called covariant by Holevo B. 

We restrict our search for global optimal measurements to POVMs that have one 
additional property: the POVM consists of múltiples of "projection operators" onto 
(generally unnormalizable) states \x), 

E(x) dx = — - \x)(x\ (35) 
(C is a real constant). The motivation for this assumption is that measurements 



not described by one-dimensional "projectors" have less resolution fï9fl , but it would 
be useful to make this motivation precise or to investigate whether covariant mea- 
surements that do not satisfy property fl35|) can be optimal. Notice that we do not 
require that the states \x) be orthogonal, and if they are not, they are necessarily 
ov ercomplete. The constant C could be absorbed into the states \x), but it is useful 
to leave it free so that these states can be given conventional normalizations in the 
examples of Section [|. 

Without loss of generality we can discard the freedom to re-phase the states \x), 
because the POVM is unaffected by re-phasing, and thus replace the displacement 
property of the POVM with a displacement requirement on the states, 

e~ ixh \x) = \x + X) . (36) 
This displacement property, written as 

(x\e- ixk \víj) = (x- X\ij)) = e~ xd/dx (x\^) , (37) 

is equivalent to saying that in the x representat ion, h is represented by a derivative: 

1 d 

h^-^r- (38) 
i ox 

The probability density fl33l) can be written as 

p(x\X) = |<Mx)| 2 = \Mx ~X)\ 2 = p(x - X) , (39) 

where 

Mx) = ~^Wx) = ^( x ~ X \^o) = M% - X) (40) 
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is the "wave function" of the state vector \ipx) in the x representation. 

Equations (|32"D , (|35|) , and ( |3~6"D are the three properties that we require of the 
POVM E(x)dx. Holevo [5(Chap. IV. 7)] considers the same sorts of measurements; 
his treatment, while more rigorous mathematically than ours, is inaccessible to many 
physicists. The three properties are preserved by a "gauge transformation," which 
replaces the states |x) with states 

e lfCh) \x) , (41) 

where / is an arbitrary real-valued function. This gauge freedom plays an important 
role, as we discuss further in Section [372] and in the examples of Section ^. 

If the POVM E(x) dx is an optimal measurement, then it saturates the second 
inequality in Eq. (B|), which simplifies to 



/ 



d * l W= F -£r^ MiAhr) - (42) 



In this inequality we put the Fisher information F in a new form, which applies to a 
covariant measurement and which is explicitly independent of X. 

For a measurement described by the one-dimensional "projectors" |x)(rr|, the nec- 
essary and sufficient condition for an optimal measurement, as shown in U, is that 

lm((^ x \x)(x\(^j^) ]=0 for all x and all X. (43) 



Using Eqs. (|2Ï|) and (p4|) and writing 



^<a#o> = M*) = r(x)e l@ ^ , (44) 



r(x) = \Mx)\ = \[p(x) , (45) 
one can recast condition (|43|) as 

= —Im{i(ip \x)(x\Ah\^ Q )) =r 2 (x)[Q'(x) - (h)} for all x, (46) 

which is equivalent to Q(x) = {h)x + constant. After discarding the irrelevant overall 
phase due to the constant, the resulting wave function is 

il>o(x) = r{x)e i{k)x . (47) 

The POVM E(x) dx thus describes a global optimal measurement if and only if the 
wave function ipo{x) of the fiducial state is (up to an overall phase) an arbitrary real 
function times a simple phase factor that accounts for the expectation value of h. For 
a fiducial state whose wave function has a phase that is nonlinear in x for all choices 
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of the states \x), we cannot rule out the existence of a global optimal measurement, 
but we can say that any measurement that satisfies properties (|32|) , fl35|), and ([36|) is 
not optimal. 

We can get at condition fl46[) directly by calculating the mean and variance of h in 
the x representat ion, again writing ipo(x) as m Eq. 



1 Q 

(h) = I dxipQ(x)-—ipo(x) = I dxp(x)Q'(x) , (4£ 



1 [dx LÍMl + / dx p(x) [e\x) - (h)f 



i dx 

<(A&) 2 )=/^|(|-i<í.))*wf i} ~ p(z) 

(49) 

This expression for ((Ah) 2 ) connects the the Cramér-Rao bound of classical estimation 
theory [first inequality in Eq. (Dl)] to the requirements of quantum theory [second 
inequality in Eq. (pj)]. A glance at Eq. ( f42|) reminds one that the first term in ((Ah) 2 ) 
is one-quarter of the Fisher information; moreover, one recognizes that for an optimal 
measurement this first term must attain its maximum value, which is the variance 
of h. Thus, for an optimal measurement, the second term in ((Ah) 2 ), which is the 
variance of Q'(x) with respect to p(x), must be zero; vanishing of the second term is 
precisely the condition (46). 



It is instructive to consider in some detail a special case of the uncertainty re- 



lation (31), because in this special case one finds the closest connection between our 
parameter-based uncertainty relations and Standard uncertainty relations. Before con- 
sidering this special case, however, it is useful to note that the mean and variance of 
the measurement result x are given by 

(x)x = J dx xp(x — X) = X + J dx xp(x) = X + (x)o , (50) 
((Ax) 2 ) x = [ dx (x - (x) x fp(x -X)= [ dx (x - (x) ) 2 p(x) = ((Ax) 2 ) . (51) 



The mean value of x with respect to the fiducial state, (x)q, globally biases the mean 
(x)x away from the parameter. The variance of x is independent of X. 

To introduce our special case, suppose that one makes iV measurements described 
by the POVM E(x) dx and that one estimates the parameter X as the sample mean 
of the data, with the global bias removed, i.e., 

1 N 

Xest = -r= J>< - <*>«>) " (52) 
iV í=l 

This estimator is unbiased, i.e., (X est )x = (%)x — (x)o = X, and thus the deviation (||) 
becomes SX = X est — X = AX cst . In addition, the efficiency of this estimator is 
independent of N, because the mean-square deviation decreases as 1/N: 

((ÓX) 2 ) = ((AX cst ) 2 ) = ((Ax) 2 ) . (53) 
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The resulting special case of the uncertainty relation (|3Ï|) is 

((Ax) 2 )((Ah) 2 } = N((5X) 2 )((Ah) 2 } > \ . (54) 

The uncertainty relation for this estimator is identical to a Standard uncertainty re- 
lation for the measurement result x, the only difference being that the states \x) are 
generally not the eigenstates of any Hermitian operator. 

Equality in the uncertainty relation requires saturating both inequalities in 
Eq. (d). Saturating the second inequality — i.e., making F = 4((A/i) 2 ) [see Eq. (|42|)]- 
means that the fiducial wave function ipo(x) has the form (|47D . Saturating the first 
inequality means that the sample mean ( |52"D is an efficient estimator for all vàlues of iV 
and, in particular, that x — (x)q, the measurement result with the global bias removed, 
is itself an efficient estimator for N = 1. We can determine the resulting conditions by 
specializing the proof of the Cramér-Rao bound to the case of a single measurement 
with x — (x)o as the estimator. We first write the mean of the estimator in the form 

Wx - W » = /<M*-MoM* = (55) 

Differentiating this expression with respect to X and using 

d 



dX 



J dxp(x — X) = — J dxp'(x — X) (56) 



leads to 

1 = — I dx (x — (x)o)p'(x — X) 



- Jdx(x- (x) - X)p'(x - X) 



dx (x — (x)o)p'{x) = — [ dxp(x)(x — (x)o) - . (57) 

J dx 

Squaring this expression and using the Schwarz inequality yields 

d ln p(x \ \ 
dxp(x)(x - (x) ) — — — j 



< 



(| dx (x - (x)o) 2 p(x)) ( / dxp(x) ( dln ^ X) ) ) • (58) 



Re-writing the expression for the Fisher information in Eq. (42) as 



p'(%)] 2 fj , s ( dlnp(x) 



2 
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shows that Eq. fl58|) is the classical N — 1 bound on the estimator x: 



((Ax) 2 ) = Jdx(x- (x) ) 2 p(x) > - . (60) 

The condition for saturating this bound, which comes from the Schwarz inequality in 
Eq. (p8[), is that 

= _ A(x _ (x)o ) p{x) oc e -A(,-(.)o) 2 /2 (61) 

where À is a constant. 

The result of these considerations is that equality in the uncertainty relation (|5 



can be achieved if and only if the fiducial wave function has the form (|7), with 
p(x) = r 2 (x) being a Gaussian. These Gaussian states are analogous to the minimum- 
uncertainty states that give equality in the Standard uncertainty relation. Thus our 
formalism of parameter-based uncertainty relations contains within itself, in the special 
case of the estimator being the sample mean, the Standard uncertainty relation and the 
associated minimum- uncertainty states. Two points deserve mention. First, for most 
generators h, there are restrictions on the form of the wave function; these restrictions, 
which are discussed in Section |3.2| and in the examples of Section |], generally prevent 
one from choosing a Gaussian for p(x) and thus mean that there are no states that 
yield equality in the uncertainty relation fl54"|). Second, the restriction to Gaussian wave 
functions is a consequence of using the sample mean as the estimator. If one allows 
other estimators, the conditions on the fiducial wave function are weaker. Specifically, 
as we have seen, in the limit of large N, where maximum-likelihood estimation is 
asymptotically efficient, the condition for saturating the uncertainty relation ( |31~D is 
that the fiducial wave function have the form (|47|) . 



3.2 The x representation 

Up till now, it has not been necessary to construct explicitly states \x) that satisfy 
the completeness and displacement properties. Such a construction depends on the 
eigenvalue spectrum of h. Suppose that we write the eigenvalue equation for h as 

h\h,a) = h\h,a) , (62) 

where we allow for the possibility of degenerades by including a degeneracy label a. 
The orthonormal eigenstates \h,a) satisfy a completeness relation 

1 = y~]\h, a)(h, a\ . (63) 

h,a 

The displacement property (|3~6p, with x = and X — x, becomes 

(h,a\x) = e~ ixh (h,a\x = 0) . (64) 
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The displacement property thus relates all the states \x) to a particular state \x — 0), 
which remains arbitrary. 

We can now ask whether it is possible to satisfy the completeness property (|32| ) 
by noting that 



í ^\x){x\=Y, I ^ e - lx{h ~ h,) l^2\h,a)(h,a\x = 0)(x = 0\h',a')(h',a'\) . (65) 
One can arrange that 



^e-^ = 5 hh> , (66) 



in which case Eq. (K5H simplifies to 

dx 

\x){x\ = 

h 



J^\ x ){x\ = lj2\ h ' a )( h ' a \ x = °)( x = \ h > a ')( h ' a> \) ■ ( 67 ) 



To make this integral equal to the unit operator requires that 

(h,a\x — 0)(x — 0\h,a') — S aa i for all h, (68) 

which can only be satisfied if the spectrum of h has no degenerades. Thus only for 
non-degenerate h can one hope to find a global optimal measurement in terms of a 
POVM described by a single real number x. An example of how to proceed for a 
degenerate h can be found in the discussion of time-energy uncertainty relations in 



Section 4.3 



We now assume explicitly that the generator h is non-degenerate, thus allowing 
us to drop the degeneracy label a from the preceding equations. The form of the 
completeness relation depends on further properties of the eigenvalue spectrum of h. 
We illustrate the procedure here for the case that the non-degenerate spectrum of 
h is discrete (nowhere dense) and that the unitary generator e~ tXh is periòdic with 
smallest period X, i.e., e~ lXh = 1 (other non-degenerate eigenvalue spectra are dealt 
with in the examples of Section ^) . This means that all the eigenvalues can be written 
as 

h = 2-KTihl 'X , rih an integer; (69) 

any discrete spectrum can be approximated in this way for X sufficiently large. The 
periodicity allows us to restrict both the parameter X and the measurement results x 
to the finite interval [-X/2,X/2). 

The completeness condition ( p2|) now becomes 

f x / 2 dx , , , , . , í' x l 2 dx 



" ^ e - ix{h - h,) 



h 
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which can be satisfied by choosing C — X and 

(h\x = 0) = é m , (71) 

where f(h) is an arbitrary real-valued function. The completeness property thus 
requires that \x — 0) have the same magnitude of overlap with all the eigenstates of 
h. 

The minimal choice, f(h) = 0, which we distinguish by underlining, leads to 
canonical states 

jiHEltyr*-*, (72) 

h 

whereas an arbitrary choice for f(h) leads to states, 

\ x ) = \h)e im e- ixh = e if{h) \x)_ , (73) 



that are a gauge transformation (PD of the canonical states \x). A gauge transforma- 
tion corresponds to the freedom to re-phase independently each of the eigenstates of 
h — i.e., to replace \h) by e l ^ h ^\h). 

The inner product of \x) and \x') is given by 

( x \ x ') = é {x - x ' )h = e 2mnh{x - x ' )/x . (74) 

h h 

These states are orthogonal — i.e., they can be given 5 function normalization with 
(x\x') = Xbix — x') — if and only if all integers are required to represent the eigenvalue 
spectrum of h; only if the states are orthogonal — i.e., all integers are present in the 
eigenvalue spectrum — are they eigenstates of a Hermitian operator. 
The x and h representations of a state are related by 

x \ij,) = ${x) = -^J2 e lxh e~ ií{h) (h\il>) , (75) 



x vx h 

e- tm (h\^} = f X/2 dxe~ ixh i){x) . (76) 

V X J —X/2 

The amplitude e~ l ^ h ^ (h\tp) is the discrete Fourier coefficient, corresponding to integer 
ri/j, of the function ip(x), which is periòdic with period X. The wave functions ip(x) are 
restricted to periòdic functions that have vanishing Fourier coefficients for the unused 
integers. By the same token, the expansion of a state in terms of the states \x), 

rX/2 d x l rX/2 

\^)= , -yl^W) = -7= / , dxij}{x)\x) , (77) 

J-X/2 X sJX J-X/2 

is not unique; one can add to ip(x) any periòdic function g(x) that has nonvanishing 
Fourier coefficients only for the unused integers, for such a function satisfies 

rX/2 

dxg(x)\x) = 0. (78) 



\pC J-X/2 
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This lack of uniqueness expresses the ov ercompleteness of the states \x). Both the 
overcompleteness and the restrictions on the wave functions ip(x) are consequences of 
the lack of orthogonality of the states \x). 

For H = 2imnl X, we can define a "displacement operator" 



D(H) = 



X/2 
-X/2 



dxe ixH E(x) 



X/2 d x 

— e lxH \x)(x\ 

-X/2 X 

Y,K„n h+ n H e^\h')(h\e-^ 

h,h> 



(79) 



which displaces eigenstates of h; i.e., 

D(H)e im \h) = e if{h+H) \h + H) 



ÍO) 



provided H is the difference in eigenvalues. Given a choice of phases for the eigenstates 
\h), the canonical states \x) are unique in that their displacement operator D_{H) 
displaces the eigenstates \h) without the inclusion of any phase factors. Notice that 
generally D(H) is not a unitary operator. For particular eigenvalue spectra of h, 
however, as in the examples of Section [|, the displacement operator acquires additional 
important properties. 

The x and h representations ( |T5|) and ( 76) of a state show that the condi- 
tion ( f£?D for a global optimal measurement, when written in the h representation, 
with h = (h) + u, becomes 



^1) 



Since the phases in the h representation can be removed by appropriate choice of the 
function f(h), this condition reduces to 



(h) + u ipo 



(h) u TpQ 



\2) 



ipo) is non-zero, 



To make this condition meaningful requires that whenever ({h) + u 

(h) — u is an eigenvalue of h. For general eigenvalue spectra of h, the condition (|8^ ) 
can be met by only a very limited class of states, since it requires symmetric excitation 
of eigenstates \h) symmetrically located about the expectation value of h. 



4 Examples of Generalized Uncertainty Relations 

We turn now to examples of generalized uncertainty relations, first dealing, in this 
section, with nonrelativistic examples and then turning, in Section [|, to Lorentz- 
invariant versions of uncertainty relations. 
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4.1 Spatial displacement and momentum 

The first example of a nonrelativistic uncertainty relation is the familiar one of spatial 
displacements X that are generated by the momentum operator p, i.e., h = p/h: 

\M = e- lX ^ h \ÍJ ) . (83) 

The uncertainty relation ([3Ï]) takes the form 

((5X) 2 M(Ap) 2 > > ^ . (84) 

Helstrom [|J and Holevo [5(Chap. VI. 2)] have presented parameter-based uncertainty 
relations for spatial displacement and momentum, and Dembo, Cover, and Thomas 
jTJ have reviewed the basis for such uncertainty relations in the properties of Fisher 
information. 

To investigate the possibilities for optimal POVMs E(x) dx = \x)(x\dx/C, start 
from the complete set of 5 function normalized eigenstates \p) of p: 

{p\p'} = 2Trh8(p-p / ) , (85) 
dp 



2wh \p)(p\- (86) 
The displacement condition fl36|), with x = and X — x, becomes 

(p\ x ) = e ~ ixp/h (p\x = 0) , (87) 

which leads to 

, ,( f°° dx . ,, ,Y ,. 2nTi\(p\x = 0)| 2 ,. . . 

(p\ ( -ç \x) (A\ W) = ^ ~S(P - P') ■ (88) 

Thus the completeness condition (|32|) can be satisfied by choosing C — 1 and 

(p\x = 0) = e if{p) , (89) 

where f(p) is an arbitrary real-valued function. In this case, because the spectrum of 
p covers the entire real line, the states 

\x) = l°° ^ \p) e if ^e- lxp/h (90) 
have 5 function normalization, 

(x\x') = 6{x - x') , (91) 
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and thus are eigenstates of the Hermitian operator 

/oo A /*oo 

dxxE(x) = / dxx\x)(x\ . (92) 
-oo J —oo 

The minimal choice, f(p) = 0, leads to the canonical position states, 

w = r ^ e ~ txp/n ' ^ 

— j-oo znri 
which are eigenstates of the canonical position operator 

x— / dxxE_(x)= / dxx\x) (x\ . (94) 



oc 



Measurements described by É_{x) are thus canonical measurements of position. An 
arbitrary choice for f(p) leads to the states (PD|) , which, written as 

\x) = e í/tf) |x) , (95) 



are seen to be a gauge transformation ( 4Ï|) of the position eigenstates. The state |x) 
is an eigenstate, with eigenvalue x, of the operator (0), which can be written as 

x = e im xe- if{p) = x + hf{p) ; (96) 

measurements described by E(x) are measurements of this operator. Notice that x 
and p satisfy the canonical commutation relation, 

[x,p] = %h , (97) 

the gauge freedom being precisely the freedom permitted by this commutator. 
The operator that displaces momentum eigenstates, 

/oo 
dxe ixP/h \x)(x\ = e i£P/h (98) 
-oo 

[cf. Eq. flTjJ], in this CcLSG cl unitary operator, acts according to 

D(pyf(p) | p ) = e if(p+n \p + P) . (99) 

The canonical states \x) lead to a displacement operator D_{P) that displaces the 
momentum eigenstates \p) without the inclusion of any phase factors. 

Writing the position wave function of the fiducial state as ipo( x ) — r(x)e l&( - x \ the 
general relations (|48| ) and(f4ÜD for the mean and variance of h become in this case 

dx^(x)-—ijj {x) = / dxp{x)hQ'{x) , (100) 

3 1 ÜX J —OO 
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<(Ap) 2 > 



dx 



h 
1 



— oo 
2 



00 , p W 

eix + 
-00 p(x) 



dx p(x)[hQ' (x) — 



(101) 



For the minimal choice [f{p) = 0] and its canonical position operator, several authors 
have drawn attention to the way the momentum variance splits into the sum of the two 
parts in Eq. (|101|) . Stam |20| noted long ago that the variance of p is bounded below 
by the Fisher information for position measurements, Cohen has discussed and 
illustrated with examples the split of the momentum variance, and Sipe and Arkani- 
Hamed have used this split and the similar split of the variance of x to contrast 
the coherence of pure and mixed states. 

The condition for a global optimal measurement is that the position wave function 
of the fiducial state have the form 

(a#o> = Vo(aO = r(x)é {í > )x/h (102) 

[cf. Eq. P?!)]. Transforming to the momentum representation, with p = (p) + u, 



im)+u) Up) + u 



dxe- ixu/h r(x) , 



one sees that the optimality condition can be written as 



-if{(p)+u). 



((p)+u\^ ) = e^- u \(p)-u\^ o y 



(103) 



(104) 



[cf. Eq. (|8Ï|)1. If one is restricted to canonical position measurements, for which f(p) 
0, the condition for optimality is that 



(p) + u Vo) = ((P) ~ u Vo) 



(105) 



If one allows gauge-transformed measurements, then the gauge transformation can 
be used to remove the phases in the momentum representation, so the condition for 
optimality is the weaker condition that 



(p) +u Vo) = ((P) ~ u ipo) 



(106) 



i.e., that the momentum probability density is symmetric about (p). 

It is instructive to illustrate these ideas with an extended example based on a 
specific fiducial state. For this purpose, introduce an "annihilation operator" 



1 í x ,Lp s 
71 \L +l ~h , 



(107) 
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where L is a constant that has dimensions of length, and a "vacuum state" |vac), 
which is the state annihilated by à, 

a| vac ) = . (108) 

One easily verifies from this equation that in the vacuum state, x and p have zero 
mean, and their covariance matrix is given by 

(vac|(A£) 2 |vac) _ L 2 (vac| (Aj3) 2 |vac) 1 

Z 2 " = r? = 2 ' (109) 

(vac|(AiA]3 + ApA£)|vac) = . (110) 

The vacuum state is thus a minimum-uncertainty state for x and p. It is convenient 
throughout the remainder of this example to use units such that L = 1, a choice that 
gives x and p/h equal variances in the vacuum state. 

The next step is to introduce the squeeze operator |23j 

S = expQr (e" 2 ^a 2 - e 2 ^ 2 )) , (111) 

which is a function of a squeeze parameter r > and a squeeze angle (p. The squeeze 
operator has the property ]23] 

SàS^ = à cosh r + a}e 2%v sinh r 

= —j= ^x (cosh r + e 2líf sinh r) + i^(coshr — e 2lLp sinh r)^j = x . (112) 

The fiducial state in this example is generated from the vacuum state by the squeeze 
operator, 

|Vo) = £|vac> , (113) 

and is sometimes called the squeezed vacuum state. An immediate consequence of the 
property (|112|) is that the squeezed vacuum state is annihilated by x: 

x\tM = . (114) 

One can get a better feel for the nature of the squeezed vacuum state and, in 
particular, its parameters r and (p by considering x and p/h to be coòrdinates on a 
phase plane and then rotating by angle (p to new canonical coòrdinates and p' /%: 

-j={x + ip/h) = à = à'e iip = -j={x! + ip'/h)é* . (115) 

In terms of the rotated coòrdinates the operator \ assumes the form 

X = e típ (a coshr + à^sinhr) = ~^ e ' líP + J ' (^6) 
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which, together with Eq. (|114j) , implies that in the squeezed vacuum state, x[ and p' 
have zero mean, and their covariance matrix is given by 

<mm')w = M^a*> e -- = i _ (117) 

{ip \(M'Ap' + Ap'Ax)\ipo} = . (118) 

The squeezed vacuum state is thus a minimum-uncertainty state for the rotated 
coòrdinates x' and p'; relative to the vacuum state, x[ has uncertainty reduced by 
a factor e~ r , and p' has uncertainty increased by a factor e r . Figure [I] depicts the 
squeezed vacuum state on a phase-plane diagram. 

If one rotates to any other orthogonal axes, the position variance gets bigger than 
the variance of x/ (recali that r > 0), because the reduced variance of x is contam- 
inated by the increased variance of p'. Indeed, the covariance matrix of the original 
canonical coòrdinates, obtained directly from Eq. ( |114| ) or by rotating back to the 
original coòrdinates, is given by |23j 

<MA£) 2 |^o> = ~ (e" 2r cosV + e 2r sinV) = , (119) 

(?/>o|(Ap) 2 |V>o) 1 / _ 2r . 2 2r 2 x 1 non x 

I^MMAp+AgA^ = _i smh2rsm 2y = _j^) = M;- 1 ) (121) 

^ 2 y 2Re( 7 ) 2Re( 7 - 1 ) ' v ; 

where 

cosh r + e 2 *^ sinh r 1 + i sinh 2r sin cosh 2r + sinh 2r cos 2(p 

7 = ; = = (122) 

cosh r — e 2líp sinh r cosh 2r — sinh 2r cos 2y9 1 — i sinh 2r sin 2ip 

is a complex constant. This covariance matrix can also be gotten from the wave 
function of the fiducial state in the canonical position representation , 

Mx) = (x\ } Po) = (^) 1/4ex p(-^) > ( 123 ) 



which follows from integrating the differential equation that represents Eq. (|114|) in 
the canonical position representation. An irrelevant phase factor is set equal to unity 
in the wave function (|123|) (Schumaker |23| has given a consistent set of phases for 
squeezed-state wave functions). 

It is now straightforward to find the optimal measurement. The wave function in 
the momentum basis is given by 

(í#o> = J^dxe-^Mx) = ^H(47rRe( 7 - 1 )) 1/4 exp('-^") , (124) 
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where J I7I/7 is an overall phase factor. According to the optimality condition ( 104 ), 
choosing f(p) to cancel the imaginary part of this complex Gaussian, i.e., 

f(p) = -h^-^/h 2 , (125) 

yields an optimal measurement, corresponding to measuring the operator 

x = x + hf'ift) = x — Im(7 _1 )p//ï, . (126) 

The distinguishing feature of using a squeezed state as the fiducial state is that the 
optimal measurement is a linear combination of x and p. Transforming to the x 
representation yields a real wave function 

7 f 1 \ 1 x • 




M*) = (*M = \l-[^^) ex p(-^fpy J > ( 127 ) 



aside from the overall phase factor \/ I7I/7, in accordance with the general condition 
for an optimal measurement. 

One feature of the optimal measurement in this case, which follows from the fact 



that ipo{x) is a Gaussian wave function of the sort considered at the end of Section |37ï 
deserves emphasis. The probability density of measurements of x, 




p(x) = \Mx)\ 2 = ( ^RlZTZn ) <M -T7-77^T ! • < 12 ' S ' 
is a zero-mean Gaussian with variance 
<MA*)U> = ÍRe( 7 -) = i (e-si„V + e-cosV)- 1 = ■ (™) 

Generally one must appeal to the large-iV asymptotic limit to saturate the first (clas- 
sical) inequality in Eq. @ — i.e., to achieve the Cramér-Rao bound — but since the 
statistics of x are Gaussian, no such appeal is necessary. Indeed, for Gaussian statis- 
tics the sample mean ( [52]) of the data (here (x)q = 0) provides an effi cient estimator 
for all vàlues of N, as is discussed at the end of Section The Gaussian statistics of 
x for the fiducial state, displaced according to Eq. (p3|), imply that (X cst ) x = X — i.e., 
the estimator is unbiased — which means that the estimate's deviation away from the 
parameter becomes SX = X est — X = AX est . The mean-square deviation is indepen- 
dent of X and reduces to 

<(5X) 2 > = <(AX est ) 2 > = l<MAx) 2 |<A)> = Jïf ■ (130) 

The final equality, which shows that X cst is an efficient estimator, follows most easily 
from the form of the Fisher information in Eq. d59|). In this case, where an efficient es- 
timator is known, one can proceed directly to equality in the uncertainty relation 1 
without going through the Fisher information, by combining Eqs. ([L29|) and (|ï"3( 
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One gains insight into the optimal measurement by writing the measured opera- 
tor pD as 

~ . , P. ü i cos 9 + (p/h) sinO 

x = x+ —ta.ia.9 = , (131) 

n cos 9 

where 

t / i\ sinh 2r sin 2g? 

tanfl = -I1117- 1 = , 132 

cosh 2r + smh 2r cos 2tp 

and regarding í as a species of position operator that arises from a rotation in the 
phase plane by angle 9, followed by rescaling by 1/ cosé*. The rescaling means that 
displacement by X produces the same "signal" in x as it does in x. The optimal 
angle 9 is not equal to tp, the rotation angle that minimizes the variance of the rotated 
position; instead, the optimal angle is a compromise between reduced "noise" and 
reduced signal, both of which come with rotation (see Fig. p]). The rescaling of x 
accounts for the reduced signal, so the variance of x, 

= (íp \(Ax) \ip ) + 2 A tanóH -3 tan 9 , 

f v n 

(133) 



is a noise-to-signal ratio p4f . Indeed, the angle that minimizes this noise-to-signal 
ratio, 

, |(^o|(AiAp + ApAi)|V;o)/^ 
(^o|(Ap) 2 |^o)/^ 
defines the optimal measurement. 

4.2 Harmonic-oscillator phase and number of quanta 

For our second example of a nonrelativistic uncertainty relation, consider a harmònic 
oscillator that has creation and annihilation operators and à. The "number opera- 
tor" 

n = à^à (135) 

has eigenstates \n), called "number states," where n — 0, 1, ... is the number of quanta. 
Shifts X = $ in the phase of the oscillator are generated by the number operator, 

|V*> = e^foM , (136) 
i.e., h = —fi. The uncertainty relation fl3T|) then reads 

(m%((Anf) > ^ . (137) 

Holevo [5(Chap. III. 9)] has considered this sort of phase uncertainty relation. Lane, 
Braunstein, and Caves [OH], in a detailed analysis of phase measurements, have used 
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the formula ([49]) , specialized to give the variance of the number operator, to bound 
the Fisher information for the phase. 

A phase shift $ can be thought of as a dimensionless time [measured in units of 
(harmonic-oscillator period) /2tt] , so the uncertainty relation (|137|) is a dimensionless 
time-energy uncertainty relation, special because of the uniform spacing of the eigen- 
states of the generator n. General time-energy uncertainty relations, corresponding to 
other energy spectra, are considered in Section |4T3} 

Since phase shifts are periòdic with period 2ir, $ can be restricted to the inter- 
val — 7r < $ < 7r. It might be thought that there is a difficulty with the phase- 



number uncertainty relation (|137 ) when the fiducial state is a number state, for which 
((An) 2 ) = 0; the uncertainty relation then forces ((óQ) 2 )q> — > oo, even though a sen- 
sible estimator $ est is restricted to the same 2ir interval as is $. No difficulty arises, 
however, because for a number state, no measurement can provide any information 
about the phase shift; thus, any estimator, sensible or not, satisfies d(§ cst )/d§ = 0, 
with the result that the deviation 5$ of Eq. flD diverges, even if <3> es t is restricted to 
a finite range. 

The possibilities for POVMs E((f>) d(p = \4>)((f>\ d(p/C (— n < <fi < tt) are determined 
by the displacement condition (ÜO), which, with = and $ = 0, becomes 



( n |0) = e i<M (n|0 = 0) , (138) 
and by the completeness condition (j32|), 

d<h ?tt 00 

1=/ 7?l0)(0l = 77El^l0 = O)| 2 H(n|, (139) 

which can be satisfied by choosing C = 2n and 

(n\<j> = 0) = e i/(n) , (140) 

where f(n) is an arbitrary real-valued function. Since there are no number states for 
negative integers, the phase states \(j>) are not orthogonal, the inner product being 
given by |26 



') = E 




i(4>—tf>')n 



^2 e -i{4>-4>> + sgnHe- í( ^' )n + 1 



^(0-0')-^ot(^)+l, (141) 



where 



sgn(n) 



•1 , n < 0, 



, n = 0, (142) 



1 , n > 0. 
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Hence the states \<p) are overcomplete and are not the eigenstates of any Hermitian 
operator. There is no Hermitian phase operator in the infinite-dimensional Hilbert 
space of a harmònic oscillator [|, |26], |27|, p9[ , although one can be constructed if 
the harmonic-oscillator Hilbert space is truncated to be finite-dimensional |30| . 

The minimal choice, f(n) =0, leads to the Susskind-Glogower |27j canonical phase 
states, 



n=0 

which are eigenstates of the non-unitary number-lowering operator 



(143) 



oí(j> 



[n 



+ 1) -l/2 



a = an 



-1/2 



53 |n- l)(n| 



71=1 



í.e., 



i4>\ 



(144) 



(145) 



Helstrom 0] and Holevo [5(Chap III. 9)] have considered measurements described by 
the Susskind-Glogower states. An arbitrary choice for f{n) leads to states, 



EH 

n=0 



D if(n) i<jm _ Jf(n)\ 



(146) 



that are a gauge transformation ([y]) of the Susskind-Glogower states. The state 
is an eigenstate, with eigenvalue e 1 ^, of the operator 



if(n) e i(p e -if(n) 



-j[/(n+l)-/(n)]^ ^ e -i[/(n)-/(n-l 



n — l)(n|e" 



n=l 



-if(n) 



(147) 



the differences f(n + 1) — f{n) and /(^) — /(ri — 1) in the exponents are analogous to 
the derivative hf'(p) in Eq. (^). 

For iV an integer the number displacement operator is given by 



D(N) 



d(f)e- i4,N È((f>) 



7T Z7T 



if(n') 



n) (n\e 



-if(n) 



]T e í/(n+Af) \n + N) (n\e- if{n) (148) 



n=max(0,-iV) 



[cf. Eq. (|79|)1. Because there are no number states for negative integers, D(N) is not 
unitary; the final form of D(N) is a consequence of the regular spacing of the number 
states. Notice that D(— 1) = e l f( n )é<t> e r % f( n ) [thus the states \<p) are eigenstates of 
D(-l)} and £>(1) = [D(-l)]ï = e if ^é^ e~ if ^ . 
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The <p and n representations of a state are related by 



(</#) = ^(0) 



(149) 



^7T 



(150) 



e -*/(n) (n|^) being the Fourier coefíicient of the periòdic function ^(0). The condition 
for a global optimal measurement, that the wave function of the fiducial state have 
the form 

M<P) = r(0)e-*W* , (151) 



is equivalent to the following requirement on the number-state amplitudes: 



</«»>+«) ^) + = e í/(<n) - u) ((n) -mU ) 



(152) 



[cf. Eq. (0)]. If one is restricted to Susskind-Glogower phase measurements [/(n) = 0], 
the condition for optimality is that 



(n) +m Vo) = ((n) -u ip ) 



(153) 



but if one allows gauge-transformed measurements, the condition for optimality be- 
comes 



[n) + u íp 



(n) -u i/j 



(154) 



In either case, the condition for optimality can only be met by a limited class of states; 
in particular, because of the discreteness of the number states, (n) must be integral or 
half-integral, and because of the lower bound at n = 0, (n\ip) must vanish for n > 2(n). 

Since the optimality conditions appear to be so restrictive, it is worth noting that 
a large class of "semiclassical" states satisfy them approximately. By a semiclassi- 
cal state, we mean one that has number amplitudes (n\ip) that are concentrated at 
large n, rendering the lower bound at n = irrelevant, and are spread over a wide 
range of vàlues of n, making the discreteness of n unimportant. For semiclassical 
states measurements described by E(4>) are nearly optimal provi ded only that the 



number probabilities 



(n 



are symmetric about (n) [cf. Eq. (|106|) 1. The extent to 



which measurements of E(<fi) are sub-optimal for semiclassical states deserves further 
investigation. 



4.3 Time and energy 

For our final example of a nonrelativistic uncertainty relation, consider the Hilbert- 
space path traced out by dynamical evolution under the Hamiltonian H: 

llM = e- iTÈ ' h \%) . (155) 
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The parameter here is the elapsed time T, and the temporal displacements are gener- 
ated by h — H/h. The uncertainty relation ( |3lD reads 



((5T) 2 > r ((A#) 2 > > ^ . (156) 

This inequality means that no matter what measurements are made to determine the 
elapsed time T and no matter how the data from those measurements are processed to 
give an estimate of T, the estimator's mean-square deviation from the actual elapsed 



time must satisfy Eq. ( 156 ) 



The time-energy uncertainty relation (|156|) must be used carefully, however. For 
example, suppose one wishes to estimate elapsed time from the dynamics of a small 
system decaying into an environment. The inequality ( |156|) places useful limits on 
such an estimate only if one uses the total Hamiltonian of the system and the environ- 
ment. An alternative approach, which focuses on the dissipative dynamics of the small 
system, is to use a master equation to describe the dynamics of the small system, to 
compute ds/dT from the master equation, and then to use the original inequality (J9[) 
to place limits on the estimation of elapsed time |JT) . 

Mandelstam and Tamm |3j derived the first parameter-based uncertainty relation, 
for time and energy, in the following way. They realized that to measure elapsed time 
T, one measures an observable A that changes with time — a clock observable. By 
defming a time uncertainty 

AT ^ ((Ai) 2 ) 1 / 2 = ft((AÍ) 2 )V 2 
\d(À)/dT\ \([À,H))\ ' 

they converted the Standard operator uncertainty relation for A and H, 

((AÀ)Y 2 ((mY">l\([ÀH}}\, (158) 



into a time-energy uncertainty relation 



AT((AH) 2 ) 1/2 > ~ . (159) 

The key idea in Mandelstam and Tamm's work, to regard elapsed time as a pa- 
rameter to be determined by measuring some other quantity, underlies the formalism 
of parameter-based uncertainty relations. The technical advances in the present for- 
malism are, first, the use of estimation theory to incorporate easily the possibility 
of múltiple measurements and to quantify precisely the precision with which a pa- 
rameter can be determined and, second, the use of POVMs to allow for all quantum 
measurements that might be used to infer the parameter. Helstrom Q and Holevo 
[5(Chaps. III. 8 and IV.7] pioneered in using these technical advances to formulate 
time-energy uncertainty relations. Hilgevoord and Uffink |2], || have formulated a 
different sort of parameter-based time-energy uncertainty relation. 
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For the case of pure-state time evolution, Anandan and Aharonov [16] noted the 
connection between the Hilbert-space angle (|27|) and the variance of the Hamiltonian 
H . This connection follows from combining Eqs. (E6f) and (f29|): 



de ids ((A#) 2 )f 

dT = 2dT = h • (160) 

Knowing that Hilbert-space angle is related to distinguishability through the inner 
product, Anandan and Aharonov formulated an uncertainty relation by asking for the 
minimum time for the evolution to proceed to an orthogonal state. Anandan WR and 



Uhlmann |32| generalized this approach to mixed states. Our formulation differs in 
that we also relate Hilbert-space angle to statistical distance and thus to a precise 
measure of the uncertainty in determining the elapsed time T, i.e., the minimum 
mean-square deviation ((5T) 2 )x- 

The states \t) that are used to describe global optimal measurements can be ob- 
tained from the energy eigenstates |e): 

H\e) = e\e). (161) 

If the spectrum of energy eigenvalues is discrete and non-degenerate, then the time 
representat ion follows immediately from obvious changes in the notation of Section || 
For example, the time states are given by 

\t) = J2 \e)e if{e) e- ite/h , (162) 

€ 

with the minimal choice, /(e) = 0, giving the canonical time representation. The 
states \t), like position eigenstates and phase states, are generally not physical states, 
as they typically have infmite energy. Holevo [5(Chaps. III. 8 and IV. 7)] has considered 
the canonical time representation and its application to optimal measurements and 
has worked out in detail the example of a free particle, where the energy spectrum 
is continuous and doubly degenerate. We review the free-particle example here to 
provide an example of how to proceed when the generator h is degenerate. 
Consider then a free particle with Hamiltonian 

H = f/2m . (163) 

The energy eigenstates coincide with the momentum eigenstates \p), which we nor- 
malize as in Eq. fl85|) . The energy eigenstates are, however, doubly degenerate (except 
for p — 0), with eigenvalues given by 

e = p 2 /2m . (164) 
A convenient way to deal with the degeneracy is to introduce a degeneracy label 

a = sgnQo) (165) 
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[\p\ = °V\ cf- Eq. (|142| )1, which allows us to write 



p = ffV2me . (166) 

The energy eigenstates can now be defined as 

\e,a) = (m/2e) 1 ' A \p) ^ \p) = {ap/mf^a) , (167) 

where a is used to distinguish degenerate energy eigenstates and where the normal- 
ization is chosen so that 

(e, cr|e', a') = 27rft<W<5(e - e') , (168) 

. . POC rjf 

1 = E/ ^\e,a)(e,a\. (169) 
a Jo 2nh 

We can now find global optimal measurements in terms of time states \t,a), where 
the states |í, +1) are constructed as in Eq. ( |162| ), but in the a = +1 subspace of 
Hilbert space, and the states \t, —1) are similarly constructed in the a — — 1 subspace. 
Notice, however, that because of the degeneracy we have the freedom not only to re- 
phase each of the energy eigenstates independently, but also to use as the bàsic energy 
eigenstates any orthonormal linear combination of the states |e, +1) and |e, — 1). In 
symbols, we have the freedom to choose new energy eigenstates 

\e, 1 ) = J2\^vy m U a ,(e), (170) 

where 7 = ±1 is a new degeneracy label and U al (e) is a 2 x 2 unitary matrix with 
unit determinant. 

With this freedom in mind, we seek a global optimal measurement in terms of a 
POVM E(t, 7) dt, where the possible results of the measurement are labeled by the 
continuous parameter t and the discrete parameter 7 = ±1. The POVM satisfies three 
properties analogous to Eqs. (|32"1), (|35l), and ([36]): 



dt 

E(t, 1 )dt = -\t, 7 )(t, 7 \, (171) 



C 



1 = W dtE{t,i)=Y, 77M>(í,7|, (172) 

^ J— 00 ^ J— 00 L·i 

e -iTÈ/n^^ = | í + Tj7 ) . (173) 
The displacement condition (|173 ), with t = and T = t, becomes 

(e, ( x|t ) 7) = e-^(e,a|í = 0,7), (174) 
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which leads to 

*l (£ % I*. 7) (t, 7i) K, ✓) = ±p(e - 6') £<e, <r|t = 0, 7 ) (t = 0, 7 |e, a'> . 

(175) 

Thus the completeness condition ( |172| ) can be satisfied by choosing C = 1 and by 
requiring that 

^(e,a|í = 0,7)(í = 0,7|e,a / ) = ^ , (176) 

7 

which, in turn, means that (e, a\t = 0,7) is a 2 x 2 unitary matrix. By removing the 
common phase factor from this unitary matrix, it can be written as 

(e,(r|í = 0,7) = e*^£/ O7 (e), (177) 

where U al (e) is the unit-determinant unitary matrix of Eq. ( |1 70| ) . Notice that the new 
energy eigenstates (|170|) satisfy (e, 7^,7') = 8 ri ie~ %u ' h . 

Because the energy spectrum is bounded below, the time states 

\t, 7) = E i" P* l e ' °y m U^{e)e-^ = r A | e , 7>e -«e/» (178) 
a Jo inri Jo zirh 

are not orthogonal, their inner product being given by 

fc*^=^f ss^^^'-^+éKrb)) ■ (179) 



The canonical time representat ion results from choosing /(e) = and U ay (e) = 8 ai . 

The probability density that a measurement yields results t and 7, given parameter 
T, is given by 

p(í, 7 |T) = \Mt,l)\ 2 = \Mt-T,l)\ 2 = p(t-T, 7 ) , (180) 

where 

Mt, 7) = (t, 7|^t> = (t — T, 7^0) = ^o(í - T, 7) (181) 

is the wave function of the state \ipT) m the time representat ion. The displacement 
property (|173[ ) implies that in the time representation, is represented by a derivative: 

* ~ "I ■ (182) 

Writing the time wave function of the fiducial state as ipo{t, 7) = r(í, 7)e ï0( ' í ' 7 \ the 
general relations (58) and (fï9 ) for the mean and variance of /i become in this case 

(#> = £/ dtr (t,lh-KiMt,j)=J2 dtp(t,>y)h&{t,>Y) , (183) 
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<(A^) 2 > = ET 



^-i(F)U (í, 7 ) 



4 ^ 



P(*,7) 



/oo 
rft^,7)[^e'(t, 7 )-(ií)] 2 . 
-00 



;i84) 



Using this expression for the variance of H or using the general condition fl43"| ) for 
an optimal measurement, one can derive that the requirement for a global optimal 
measurement is that 0(í, 7) = (H)t/h + (constant)^. Discarding an irrelevant overall 
phase due to the two constants, but retaining the differential phase between the 7 = ±1 
parts of the wave function, one can write the resulting fiducial wave function for a 
global optimal measurement as 



^(í,7) = e íw Kí>7)e í<A>í/ ", 

where \i is a constant. 

The time and energy representations of a state are related by 



(185) 



(í,7l 



^7) =E/ ° O ^/V^)[/; 7 (e)(e,a| 



00 de 



2nh 



e,7 



/OO 
díe- ifc /V(í,7)- 
-00 



(186) 
(187) 



These representations show that the condition (|185|) for a global optimal measurement, 
when written in the energy representation, with e = (H) + u, becomes 



(H) +M,7 7p 



Y,e~ lfi{È)+u) U;J(H) +u){(H) +u,a ^ 



E 



Jf((H) 



-^U ai {{H)-u){{H)-u,a\^ 



(H) -u,7 ^0 



where we discard the differential phase /1 because it can be absorbed into the unitary 
matrix U ai . The condition ( |188|) can be satisfied, by appropriate choices for the 
function /(e) and the unitary matrix U aj (e), if and only if 



(189) 



i.e., the total probability density to have energy (H) + u is the same as the total 
probability density to have energy (H) — u. 
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5 Lorentz-Invariant Uncertainty Relations 



We now apply the formalism developed in Section [2.2| to formulating Lorentz-invariant 
uncertainty relations for estimation of the displacement and Lorentz-rotation param- 
eters of the Poincaré group. We deal first with the displacement parameters, where 
we are seeking a restriction on the estimation of a space-time translation and, hence, 
on the estimation of the invariant space-time interval. The generator of space-time 
translations is the operator for the energy-momentum 4- vector 

p = p a e a = P°e + P = P°e + P j ej , ( 190) 

for whatever fields are used to distinguish translated frames. We write the displace- 
ment 4-vector as 

X = Sn = Sn a e a , (191) 

where 

n = n°e + n = n°e + n 3 ej (192) 

is a (timelike or spacelike) unit 4-vector that gives the direction of the space-time 
translation and S is the invariant interval that parametrizes the translation. The path 
on Hilbert space is given by 

\il>a) = e~ íSn - p/h \^ ) , (193) 

where 

n P = T] aP n a P p = n a P a = -n°P° + n-P . (194) 
Here ||?7a/3|| — diag(— 1, +1, +1, +1) is the Minkowski mètric of special relativity (we 

adopt units such that the speed of light c = 1), and n ■ P — n^P^ is the three- 
dimensional dot product. 

From Eq. ([H]) the uncertainty relation for estimation of the invariant interval S is 

((^) 2 ) 5 ((n • AP) 2 ) = ((5S) 2 ) s n^(AP a AP p ) > ^ . (195) 

When n is timelike, this is a time-energy uncertainty relation for the observer whose 4- 
velocity is n, and when n is spacelike, this is a position-momentum uncertainty relation 
for an observer whose 4-velocity is orthogonal to n. In particular, when n = eo, the 
time-energy uncertainty relation takes the form 

<(^) 2 M(AP°) 2 > > ^ , (196) 

and when n = n = n^ej is a spatial unit vector, the position-momentum uncertainty 
relation becomes 

((6Sr) s ((n-Apf}>^. (197) 
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For illustration, suppose that the relevant field is the free electromagnètic field. 
When considering the energy-momentum 4- vector as a generator, it is most convenient 
to decompose the field in terms of plane-wave field modes, for then 

p = E^ k 4,A> ( 198 ) 

k,a 

is a sum of separate contributions from the various modes. In Eq. ( |198| ) k = ue + k = 
ueQ + Wej is a (null) wave 4- vector, with u = \k\ = k; the sum runs over all plane- 
wave field modes, i.e., over all all wave 3- vectors k and over the two helicities, denoted 
by o. Since the generator n • P/h for any space-time translation is determined by 
the number operators for the plane-wave field modes, global optimal measurements 
will involve measurements of phase shifts of these modes. This is not a surprising 
conclusion because the effect of a space-time translation is to shift the phase of each 
plane-wave field mode. Indeed, if only a single plane-wave field mode is excited, the 
discussion of global optimal measurements of the invariant interval reduces to the 
analysis of phase measurement in Section EO. If many modes are excited, as in a 



pulse of electromagnètic radiation, the discussion of global optimal measurements is 
more complicated. Measurements of phase shifts in the multi-mode case are only 
beginning to be considered |33|, ^4[ |35fl . Notice that when many modes are excited, 
the generator n-P/h becomes highly degenerate, a situation that cannot be addressed 
by the general considerations of Section |3|. 

Turn now to the case of Lorentz transformations, where we seek restrictions on the 
estimation of the parameters that describe boosts and spatial rotations. The generator 
of Lorentz transformations is the operator for the antisymmetric angular-momentum 
two-tensor 

J = J Q/3 è Q <g> èp , (199) 
whose components are given in terms of the stress-energy tensor by 

j 1 * = J d 3 x (x a f^ - x p f a0 ) . (200) 

The path on Hilbert space is given by 

|Ve> = exp(-^eO^J^) |Vo> , (201) 

where G is the Lorentz-rotation parameter and fl al3 is an antisymmetric two-tensor 
that specifies the sense of the Lorentz rotation. 

For a boost with velocity v in the direction of a spatial unit vector n = n^e,-, B 
is the velocity parameter corresponding to v, i.e., cosh© = (1 — f 2 ) -1 / 2 , and the only 
non-zero components of Í2 a/3 are the time-space components Í2 - 7 = — fP = — nP . The 
path on Hilbert space becomes 

\ipe) = e ï0H ^ h \^ o ) , (202) 
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where the boost generator, 

K = K j ej = J j ej , (203) 

is an energy-weighted position operator. For a spatial rotation about the spatial unit 
vector n = n^ej, is the rotation angle, and the only non-zero components of íl a/3 
are the space-space components íl^ = e^ kl n l , where é kl is the three-dimensional Levi- 
Civita tensor. The path on Hilbert space becomes 

|V>e> = e-™^^) , (204) 
where the generator of spatial rotations, 

f=J>e j = y kl j kl e j , (205) 

is the angular-momentum operator. 

The general form of the uncertainty relation for estimation of the Lorentz-rotation 
parameter B is 

((50) 2 )e \n a ^(AJ af3 AJ^) > . (206) 
For a boost the uncertainty relation, 

((50) 2 )e((n-Aie) 2 )>^, (207) 

expresses the quantum-mechanical limitations on determining the velocity parameter 
O. This uncertainty relation is complementary to the relativistic position-momentum 
uncertainty relation ( |1 97| ) . In Eq. ( |197| ) the parameter is a spatial displacement, and 
the operator is the component of 3-momentum which generates the displacement. 
In Eq. ( |207| ) the parameter is related to a velocity change, and the operator is the 
component of energy-weighted position which generates the change in velocity. For a 
spatial rotation the uncertainty relation, 

((50) 2 )e((n·AÍ') 2 )>^, (208) 

expresses the quantum-mechanical limitations on determining a rotation. 

To investigate global optimal measurements of a spatial rotation or a boost, it 
would be wise to decompose the relevant field in terms of angular-momentum modes 
or "boost modes." Such an investigation lies outside the scope of the present paper. 



6 Conclusion 

Much ink has been devoted to the problem that many quantities of physical interest, 
such as time or harmonic-oscillator phase, though determined routinely from mea- 
surements, cannot be accommodated within the conventional quantum-mechanical 
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description of measurements, because such quantities have no associated Hermitian 
operator. The aim of this paper is to show that this problem is only apparent. We 
eschew tedious discussions of the status of such quantities in quantum theory. Instead 
we develop a formalism that allows us to derive quantum-mechanical limitations on 
the determination of such a quantity, without ever having to introduce an operator 
associated with the quantity, and we illustrate the formalism with numerous examples. 

The formalism is founded on the idea that such a quantity should be treated 
as a parameter, to be determined from the results of measurements. To derive strict 
quantum-mechanical limits on such a determination, we must be able, first, to describe 
all measurements permitted by the rules of quantum mechanics — this is accomplished 
by using the formalism of POVMs — and, second, to set bounds on all possible ways of 
estimating the parameter from the results of the measurements — this is accomplished 
by appealing to the Cramér-Rao bound of classical parameter-estimation theory. The 
resulting quantum-mechanical limitations are expressed as Mandelstam-Tamm uncer- 
tainty relations involving the precision of the parameter estimation and variance of the 
operator that generates changes in the parameter. These uncertainty relations take 
into account naturally the expected improvement in determining the parameter as 
one is allowed to make measurements on an increasing number of identically prepared 
systems. Moreover, we are able to derive general conditions for optimal measurements 
that can achieve the lower bound in the uncertainty relation, although it is generally 
not known how to perform such optimal measurements. The final result is a formalism 
that increases considerably the scope and power of uncertainty relations in quantum 
theory. 
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Figure 1: Phase-plane representation of optimal measurements of displacement of a 
squeezed vacuum state. The squeezed vacuum state — S\vac) is represented by 
a sòlid "uncertainty ellipse" centered at the origin. The principal axes of the ellipse 
are oriented along the directions defined by the uncorrelated coòrdinates x' and p', 
which are rotated by an angle ip relative to the axes defined by the canonical position 
x and the momentum p; the principal radii of the ellipse are given by the uncer- 
tainties <^ 1 ( Ax') 2 |^o> 1/2 = and {^(kp') 2 ^) 1 ' 2 /% = e r /V2. The dotted 

uncertainty ellipse depicts the state obtained by displacing the squeezed vacuum state 
a distance X along the x axis. The global optimal measurement for distinguishing 
displaced squeezed states corresponds to measuring a variable x [see Eq. ( |131|) 1, which 
is a rescaled position variable along an axis rotated by angle 9 from the axis of the 
canonical position variable x. The optimal measurement represents a compromise be- 
tween maximal "signal," which would be obtained by measuring the canonical position 
variable x, and minimal "noise," which would be obtained by measuring the rotated 
position variable xf [see Eq. ( |134| )]. 
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FIG. 1. of "Generalized uncertainty relations: Theory, examples, and Lorentz invariance," 
by S. L. Braunstein, Carlton M. Caves, and G. J. Milburn, submitted to Ann. Phys. 
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